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We have analyzed numerically the response and noise- 
limited charge sensitivity of a radio-frequency single-electron 
transistor (RF-SET) in a non- superconducting state using the 
orthodox theory. In particular, we have studied the perfor- 
mance dependence on the quality factor Q of the tank circuit 
for Q both below and above the value corresponding to the 
impedance matching between the coaxial cable and SET. 



I. INTRODUCTION 

An important drawback of the conventional single- 
electron transistor^'^ (SET) is its relatively large out- 
put resistance which should be much larger than the 
quantum resistance Rq = h/Ae^ ~ 6.5 krj. This lim- 
its the operation frequency of the prospective integrated 
single-electron circuits'^''* and imposes a severe frequency 
limitation for individual SETs used nowadays as elec- 
trometers. Estimating the total capacitance of a wire 
delivering SET signal from a cryostat to outside elec- 
tronics as 1 nF, we get the time constant on the order 
of 1 nF X 10'' 51 = 10^** s; therefore the operating fre- 
quency is limited to few kHz, which is a typical value 
achievable by conventional SET setups.^ While the oper- 
ating frequency can be significantly increased by placing 
a preamplifier in a close vicinity of the SET,^''' the pre- 
ferred at present solution of the problem is the use of 
the radio-frequency SET* (RF-SET) which in many in- 
stances has already replaced the traditional SET setup. 

The principle of the RF-SET operation is somewhat 
similar to the operation of the radio-frequency supercon- 
ducting quantum interference device^ (RF-SQUID) and 
is based on the microwave reflection*'^°~^^'^^'^^'^*'^^ from 
a tank (LC) circuit containing the SET (Fig. 1); another 
possibility is to use the transmitted wave.^°'^^ The mea- 
sured SET input signal changes effective SET resistance 
and affects the intensity and the phase of the reflected 
(or transmitted) wave that is later sensed by either ho- 
modyne detection or simple rectification (to separate in- 
coming and reflected waves a directional coupler can be 
used). The high operation frequency of the RF-SET is 
due to the signal propagation by the microwave, so that 
the SET does not need to charge the whole output wire, 
while the tank circuit tuned in resonance works as an 
impedance transformer providing a better match between 
the effective SET differential resistance Rd 10^ f2) and 
the microwave cable wave impedance Rq (typically 50 fi). 

For a good matching the "unloaded" quality factor 
Q = \J [Lt I Ct) I Rq of the tank circuit consisting of in- 
ductance Lt and capacitance Ct should be comparable 



to ^ Rd/Rvi ^ 50. While a much lower value, Q — 
was used in the first experiment,® the values close to 
the matching condition (sometimes even higher) are typ- 
ically used at present. Increase of the Q-f actor obviously 
decreases the RF-SET bandwidth limited by ~ u}/2Ql 
where w « 1/ \JLtCt is the "carrier" microwave fre- 
quency and Ql is the "loaded" Q- factor, which also takes 
into account the effect of the SET (see below). The 
straightforward design® shown in Fig. 1 can be somewhat 
modified to reduce the bandwidth further for a dense 
multiplexing;-'^'* however, in this paper we will consider 
only the original design. 

The RF-SET bandwidth as wide as 100 MHz has been 
demonstrated® using a relatively high carrier frequency 
ujII-k = 1.7 GHz and relatively low Q-factor Q = 6. 
However, to improve experimental RF-SET sensitivity, 
it happens to be advantageous to reduce the carrier fre- 
quency by few times (to reduce the noise contribution 
from the amplifier) and also increase Q (closer to the 
impedance matching regime), so that a practical band- 
width at present is about 10 MHz (for example, the band- 
width of 7 MHz for the carrier frequency of 332 MHz has 
been reported in Rcf.*^). 

The high operation frequency of the RF-SET makes 
it easily possible to avoid the 1// noise limitation of 
the SET sensitivity which is typically dominant at fre- 
quencies / ^ 10^ Hz, and so work in the region of 
the shot noise limited sensitivity. ^^'^^ Though experi- 
mentally the contribution from the amplifier noise is 
still comparable or larger than the SET noise, a rela- 
tively rapid improvement of the RF-SET charge sensi- 
tivity from 1.2 X 10"^ e/\/Hz at 1.1 MHz in the first 



experiment® to the value 3.2 x 10~°e/vHz (4.8 h in 
energy units) at 2 MHz reported in Ref.*^ assures that 
the pure shot noise limit will be achieved pretty soon. 
(The above numbers can be compared with the sensi- 
tivity 2 X 10~^e/\/Hz at 4.4 kHz of a purely conven- 
tional SET reported in Ref.^, the value 8 x 10"^e/\/Hz 
at 10 Hz for the "stacked" SET,24 and the sensitivity 
6 X 10-^e/%/lfe at 45 Hz for the SET made of carbon 
nanotubes.^^) 

One of important potential applications of the RF- 
SET is for the readout of the charge qubits in a solid- 
state quantum computer. **'*^~*^'^^ The possibility of a 
single-shot qubit measurement requires fast enough dis- 
tinguishing between two charge states to avoid significant 
qubit evolution during measurement. This requires suf- 
ficiently wide RF-SET bandwidth and most importantly 
good enough charge sensitivity. The estimates*^ show 
that the single-shot qubit measurement is almost within 
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the reach of present-day RF-SET performance; however, 
rehable measurement still requires significant improve- 
ment in sensitivity. This makes very important the ques- 
tion of ultimate (theoretical) RF-SET sensitivity. The 
ultimate RF-SET sensitivity is also a crucial parame- 
ter for monitoring quantum dynamics of nanomechanical 
resonators^''"^^ and for a variety of RF-SET applications 
as an electrometer in classical single-electron devices. 

In spite of significant experimental activity on RF- 
SETs, we are aware of only few theoretical papers on 
the RF-SETs. The basic theory of the shot noise limited 
charge sensitivity of the RF-SET has been developed in 
Ref.^". A similar theory has been applied to the sen- 
sitivity analysis for the RF-SET-based micromechanical 
displacement detector. ^^^^'^ Some theoretical analysis of 
the transmission-type RF-SET can be found in Ref.^^. 
The theory of a somewhat related device, radio- frequency 
Bloch-transistor, has been developed in Ref.^^. In our 
opinion, the RF-SET definitely requires further theoret- 
ical attention, since many questions about RF-SET per- 
formance has not yet been answered theoretically. 

In this paper we extend the theory of Ref.^° to the 
case of arbitrary Q- factor of the tank circuit, removing 
the assumption (strongly violated in the present-day ex- 
periments) of Q being much smaller than the impedance 
matching value. We calculate the response and sensi- 
tivity of the normal-metal RF-SET and optimize these 
magnitudes numerically over the rf wave amplitude and 
the SET background charge. Then we study the depen- 
dence of the optimized RF-SET response and sensitivity 
on the tank Q-factor, operation temperature, SET re- 
sistance, carrier frequency, and SET asymmetry due to 
asymmetric biasing. Some results of this paper have been 
presented earlier in a short form.^^ 



II. MODEL AND CALCULATION METHODS 

The schematic of the RF-SET used in our analysis is 
shown in Fig. 1. The SET consists of two timnel jimctions 
with capacitances Cij and C2j and resistances Ri and 
R2. The SET is coupled via gate capacitance Cg to the 
measured charge source (for example, a single-electron 
box or a similar structure) , which is characterized by the 
charge qs and capacitances Csi and C52 to the SET 
leads (the total source capacitance is Cs = Csi + Cs2)- 
Assuming constant qs (neglecting backaction from the 
SET), it is easy to show that the SET coupled to the 
charge source is equivalent to the simple double-junction 
SET structure with parameters 
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FIG. 1. Schematic of the RF-SET. The current I{t) 
through the SET (two tunnel junctions with capacitances C\j 
and C2i) affects the quality factor of the tank circuit (Lt and 
Ct) and therefore the amplitude and phase of the reflected rf 
wave propagating along coaxial cable (-Ro). The change of the 
source charge qs changes the eff'ective SET background charge 
go and can be monitored via measurement of the reflected rf 
wave. 



where Ci and C2 are effective junction capacitances (the 
total SET island capacitance is C^. = Ci -I- C2), and the 
total induced charge of the SET island go is the sum of 
the initial background charge (700 and the contribution 
from the measured charge source. For numerical results 
we have used the "orthodox" model^ for a normal-metal 
SET (see Appendix). 



Linear analysis of the reflected wave 

The current I{t) through the SET affects the qual- 
ity factor of the tank circuit consisting of inductance 
Lt and capacitance Ct, while the contribution to the 
tank circuit capacitance is neglected, assuming Ct ^ 
max[Cij. C2j, C5'iCs2/(Csi -I- Cs2)\- (This condition is 
usually well satisfied experimentally; it also allows us to 
neglect the effect of single electron jumps on the tank cir- 
cuit oscillations.) The quality factor of the tank circuit 
is also afi'ected by the wave impedance Rq of the cable. 
For a simple linear analysis let us substitute the SET by 
an effective differential resistance Rd- Then, assuming 
the case of weak damping, we can simply add the con- 
tributions from i?o and Rd, so that the total ("loaded") 
quality factor of the tank circuit is 



Ql = (i/Q + i/Qset)-\ 

Ro Lt I Ct 



(4) 
(5) 



where the "unloaded" quality factor Q corresponds to 
absence of the SET, Rd = 00, while Qset corresponds 
to the damping by the SET only, _Ro = 0. Notice that Q 
is fixed by the RF-SET design, while Qset and therefore 
Ql depend on the operating conditions; so even though 
Ql is a more physically meaningful quantity than Q, in 
this paper we consider unloaded Q as an independent 
parameter and call it as Q-factor. 

For the incoming voltage wave exp(icijf), the re- 
flected wave aVin exp{iujt) is determined by the complex 
reflection coefficient a: 
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Z — Ro 
Z + Ro 



Z = iwLn 



1 



iujCt + 1 /Rd ' 



(6) 



Since the measured charge signal changes the SET re- 
sistance Rd, the RF-SET response is proportional to 
d{aVin) / dRd- Using the first order approximation close 
to the resonant frequency ujq — I/^Lj-Ct of the tank 
circuit: 



Rd Wo 



(7) 



(in particular, this approximation neglects the shift of 
the resonant frequency, which is a second-order effect) it 
is easy to get 



da 
dRd 



-2Ro 



2 ■ 



(8) 



For fixed values of Rq and Rd, this expression increases 
with Q for small Q, thus showing the need of Q-factor 
increase to achieve a better RF-SET response. The max- 
imum is reached at 



Q = \J Rd/ Ra , 



(9) 



which is the case of practically matched impedances, 
Z Ki Ro,at frequencies close to resonance, and also corre- 
sponds to the condition Q = Qset = ^Ql- Notice that 
in this case the reflection practically vanishes, a ~ 0. 

Actually, Eq. (8) is not really relevant, because the rf 
amplitude Vin should depend on Q in order to maintain 
approximately constant amplitude |Vt| of the SET bias 
voltage oscillations, which is determined by the voltage 
scale of the Coulomb blockade. Taking into account the 
relation 

Vin/Vb = {Z + Ro){itoCT + l/Rd)/2, (10) 
we obtain 



da 
dRd 



-Rq Q |l + 2iQLAa>/cjo| 
l+K&(l + 2ig,Ac./c.o)^ 



(11) 



which is somewhat similar to Eq. (8) and also reaches 
maximum at Q = \/Rd/RQ- Notice that this condition 
optimizes only the RF-SET response, while the shot- 
noise-limited sensitivity can (and as will be seen later 
does) have completely different dependence on Q. 



Full analysis of the reflected wave 

The linear analysis discussed above can be used only 
as an estimate because of the significant nonlinearity of 
the SET current- voltage (7 — V) dependence. For a more 
exact analysis"^" wc use the Kirchhoff's rules taking into 
account the current I{t) through the SET. Let us sep- 
arate the voltage Vq + v{t) at the end of the rf cable 



into the dc component Vq (which can be supplied via 
the bias-tee*) and rf component v{t) = Vin{t) + Vout{t), 
where Vin{t) = Vin cosut is the incoming wave (from now 
on we do not use complex representation) and Vout{t) is 
the outgoing reflected wave. The differential equation for 
the rf component v{t) is 



w/Wq -|- v/QuJq + V = 2(1 — (J^ /wQ)Vin COSljjt 



Mm - {1)1 



(12) 



where (/) is the current through the SET averaged over 

time much longer than UJ~'^ , and the time dependence of 
the SET current I{t) can be found self-consistently from 
the time dependence of its bias voltage 

Vb{t) = Vo+ v{t) + [2Viniosm.ujt + v{t)\ Q/ujo- (13) 

In this paper we assume that the rf frequency is small 
compared to the frequency of electron tunneling through 
the SET, oj -C I/e, so that the SET shot noise is a small 
contribution compared to the deterministic part of the 
SET current I(t), which is calculated using dc I — V 
curve and Vb{t) (actually, it is still OK if this condition 
is not satisfied during some fraction of the period due 
to Coulomb blockade, since the small current does not 
affect the oscillations significantly). 

In a steady state (assuming that the induced SET 
charge qo does not change with time) the reflected wave 
can only contain the incoming frequency lu and its over- 
tones: 

oo 

Vout {t) = — Vin COS ut + [Xn COS nojt + Yn sin nujt] . 



n=l 



(14) 



We separate the term —Vin cos ut mainly to follow the 
notations of Ref.'^", even though it really makes sense 
only for a low-Q case at a; « wqj when all Xn and Yn are 
small. Using the substitution v{t) = X^^i[-^7i cos nwt + 
l^sinno;*] in Eq. (12), we find the coeflacients X„ and 
Yn as 



— RqQ 



nu) an — (3(1 — n w ) 6„ 



n2w2 + Q2(l - n'^Co^f 



Yn 



-RoQ 



nui bn + Q(l — n^cl)^) a„ 



+- 



n2(I;2-^Q2(i_„2^2)2 
2Qd;(l - 



(15) 



(16) 



where 



2 {I (t) sin nut), bn = 2 {I (t) cos nujt) (17) 



(averaging is over the oscillation period), oj = (jj/ojq is 
the normalized frequency. Sin is the Kronecker symbol. 
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and the current I{t) is calculated self-consistently using 
the SET bias voltage 



Vb{t) = Vo + 2Q — Vin sinwi + V 

n— 1 



Qnuj . . 
X cos nut + [Yn ^n) smuLot 

Wo 



(18) 



Because of the resonant behavior of the tank circuit 

at Q S> 1, the contribution of overtones (n > 2) in the 
reflected signal is always small if a; w (since the over- 
tones arc far from resonance). This can be easily seen 
from Eqs. (15)-(16), especially in the case w = ojq, when 
they are significantly simplified: 



— RqQ 



Yn — —RqQ 



nan - Q(f - n'^)K 
n2 + (52(i_n2)2 

nbn + (3(1 — rP')a„ 



(19) 

(20) 



Therefore, a linear (one-frequency) approximation in 
which only Xi and Yi are taken into account, works very 
well in this case. In our numerical analysis we used the 
result of the linear approximation as a starting point of 
the iterative procedure [V{,(^) (onj^n) {XmYn) — > 
Vh{t) —>...] to solve selfconsistently Eqs. (15)-(18) tak- 
ing into account few (typically 3-5) overtones. We have 
checked numerically that the account of overtones typi- 
cally gives a small correction in the case of a reasonably 
large Q-f actor and w « wq- 

In the linear approximation the SET bias voltage has 
only one frequency component: Vb{i) = Vo + A;, sin(cjt -|- 
(/>) and therefore for the calculation of ai and bi [see Eq. 
(17)] the SET can be simply replaced by the effective 
resistance 



Rd 



J^^ I{Vo + Ab sin x) sin x dx 



(21) 



where I(V) is the SET current-voltage dependence. [No- 
tice that /g ^ I{Vo + Ab sin x) cos x dx = 0, i.e. there is no 
effective reactance contribution.] 

Hence, this linear (one-frequency) approximation is 
completely equivalent to the case of a resistor instead 
of the SET, considered in the previous subsection. The 
only new condition is a selfconsistent relation between 
the effective resistance Rd and the amplitude Ab of the 
SET bias voltage. The amplitude Ab (which depends 
on Rd) can be calculated either using Eq. (10) in which 
^6 = 1^1 or using Eq. (18), which gives 

Ab = v/[0t:^(2y.„ - Xi) + ri]2 + {Xi + QCoYif, (22) 

while the components Xi and Yi are given by equations 



2Q^ ^- + (l-^-)^^±g^ 
Q'^ + R ^2 , Q^[i+{i-a'^)W 



(23) 



2^Q - '3^f^ff-f)^] +(i-^2)^ 

^ _ Q+R y ^24) 

where R= Rd/Ro ^ 1 (these equations can obviously 
be rewritten in a shorter way; however, they become less 
transparent to analyze). 

The equations (23)-(24) significantly simplify in the 
case (jO = wq: 



2Q\Q^ + R) 

^1 — T^TT; ^TT; ~ ^in, 



(02 + R)2 + Q2 
(Q2 + R)2 + Q2 



(25) 
(26) 



from which it is clear that |Yl/Xi| <C 1 iov Q ^ 1 , and 
therefore 

Vb{t) « Vb + sinwt = Vq + 2QLVin sinut. (27) 

Q-' + R 



Response and noise-limited sensitivity 

So far we have implicitly assumed that the SET cur- 
rent I{t) depends on time only because of the periodic 
time dependence of the SET bias voltage Vb{t). However, 
I(t) has also a small noise component, the magnitude of 
which depends on the bias voltage and therefore also has 
a periodic time dependence. The shot noise of the SET 
current leads to the fluctuations of the parameters a„ 
and bn defined by Eq. (17) and consequently to the fluc- 
tuations of the reflected wave quadratures X„ and Yn. 
Since the noise of X„ and Yn can be meaningfully dis- 
cussed only at frequencies less than lo/Ql, which is much 
less than the typical frequency 7/e of electron tunneling 
in the SET, it is sufficient to consider the low- frequency 
limit of the SET shot noise. 

The low- frequency spectral densities of a„ and 6„, and 
their mutual spectral density can be calculated as 



San = 4(5/(i) sin^ nujt), Sbn 

San.bn = 2{Sl{t) SV[\2nujt) , 



4:{Si{t)cos^nojt), (28) 
(29) 

where the averaging is over the oscillation period, Si{t) 
is the low- (zero-) frequency spectral density of the SET 
current (see Appendix), and the time dependence comes 
from oscillating bias voltage Vb{t). Consequently, the 
low-frequency spectral densities of Xn and Y^ fiuctua- 
tions, and their mutual spectral density are 

Sxn = c\{Si{t) siv? nujt) + dl{Si{t) cos^ nut) 

—Cndn {Si {t) sin 2nijjt) , (30) 

Svn = dl{Si{t) sin^ nut) + (^{Si{t) cos^ nut) 

+Cndn{Si{t) sin 2nut), (31) 

Sxn,Yn = Cndn{Sl{t) C0s2nut) 

1 



+^{dl - cl){Si{t)sm2nut), 



(32) 



4 



where 



dn 



2RoQnu} 



2RoQ^{l-n^u^) 



(33) 
(34) 



Notice much simpler equations for Sxi and Syi in the 
case CO = uq: 



Sx, =4RlQ^Si{t)sm^L0t), 
Sy^ = ARlQ'^{Si{t) cos^ ut). 



(35) 
(36) 



The change of the measured charge qs can in princi- 
ple be monitored via the change of any quadrature X„ 
or Yn. The corresponding responses can be defined as 
the derivatives dXn/dqs or dYn/dqs, while the corre- 
sponding sensitivities (minimal detectable charge 6qs for 
a small measurement bandwidth A/) are 



5qs,xn = 



\dXn/dqs\ ' 



Sqs,Yn = 



\dYn/dqs\ 



(37) 



Numerical calculations show that in the usual case 
w « Wq thus defined RF-SET sensitivity for overtones 
(n > 2) can be comparable to the sensitivity using the 
carrier frequency (n — 1). However, because of rela- 
tively small amplitude of reflected overtones in this case, 
their monitoring is impractical, and so we mainly con- 
sider monitoring of Xi and 11 which are referred below 
as X and Y. 

Monitoring both quadratures X and Y simultaneously, 

one can improve the sensitivity compared with monitor- 
ing of only one quadrature. It is easy to show that the 
resulting sensitivity can be obtained as the optimization 
over angle ip of the sensitivity corresponding to monitor- 
ing the linear combination X* — X cos ip+Y sin (p (exper- 
imentally, this is just a phase shift in the homodyne de- 
tector; notice that the contribution —Vin cos ut is noise- 
less by assumption) . The optimum sensitivity is achieved 
at tamp = [SxdY/dqs - SxYdX/dqs)/{SYdX/dqs - 
SxYdY/dqs) and the resulting sensitivity is 



X' 



SxSy — SxY 



SY{j^r+sx{£-r-2s. 



XY 



dX dV 
dqs dqs 



that can be rewritten as^° 



5qs,X' 



1/2 



(38) 



1/2 



{5qs,x) ^ + {5qs,Y) ^ - 2K/6qs,x5qs,Y 



(39) 



where K = Sxy/\/SxSy sign[{dX/dqs){dY/dqs)] is the 
noise correlation factor [Sxy is a real magnitude be- 
cause we consider only low-frequency sensitivity; for 
finite frequency Sxy in Eq. (38) should be replaced 



by KeSxy]- Notice that the response optimization 
for X* monitoring is achieved at different phase shift; 
ta,np = {dY / dqs) / {dX / dqs) that leads to the response 
dX*/dqs = [{dX/dqs)^ + {dY/dqs)^Y'\ 

Let us also consider the case when the reflected wave is 
monitored by simple rectification. Assuming for simplic- 
ity the monitoring of only the first harmonic amplitude 
Ai = [{X - Fi„)2 + y2]i/2 (overtones are filtered out), 
the RF-SET sensitivity can be calculated as 



^qs, 



[{X - Vir,)^Sx + Y^Sy + 2{X - Vin)YSxYY'^ 



\{X-Vir,)dX/dqs + YdY/dqs\ 



(40) 



while the response is obviously dAi/dqs — \(X — 
Vin)dX/dqs + YdY/dqs]/Ai. The formulas for moni- 
toring of nth overtone are similar, except Vin does not 
contribute to the amplitude An- 

In the most interesting for practice case lu = ujo the 
magnitude of Y is small in comparison with X for Q ^ 1 
[see Eqs. (25)-(26)], and the correlation factor K also 
vanishes [because di = and the SET bias voltage con- 
tains mostly sin-component - see Eq. (27)]. Then Eq. 
(39) for the best homodyne detection practically coin- 
cides with Eq. (40) for rectification and reduces to the 
sensitivity Sqs,x / V^f = V'^x /{dX/dqs) for monitoring 
of X component only (similarly, the formulas for response 
also practically coincide). Because of that, the numerical 
results for RF-SET response and sensitivity in the case 
u} = luq will assume monitoring of X quadrature. 

Since the increments of the measured charge qs and the 
effective SET charge qo are related by a constant factor. 



Sqs = Sqo{l + Cs/Cg), 



(41) 



the RF-SET response and sensitivity in respect to qs and 

in respect to qo differ by the factor 1 + Cs / Cg . All our 
numerical results will be in terms of qo measurement, so 
for X monitoring we will use the derivative dX/dqo as a 
measure of the RF-SET response and the magnitude 



5qo 



\dX/dqo\ 
as a measure of the sensitivity. 



(42) 



III. NUMERICAL RESULTS 

We have studied numerically the dependence of the 
RF-SET response and sensitivity on various parameters, 
which include "fixed" parameters (which cannot be eas- 
ily changed in an experiment) and the parameters of 
the operation point. The fixed parameters are: effective 
junction capacitances of the SET Ci and C2 (we assume 
Ci = C2 unless mentioned otherwise), resistances i?i and 
i?2 (we always assume Ri = R2), temperature T, cable 
wave impedance Rq (we always assume i?o = 50 f2), the 
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FIG. 2. Contour plots of (a, c) the RF-SET response 
dX/dqo (in units C^^) and (b, d) the noise- limited sensitivity 
5qo /\/A7 (in units ey/R-sCs) on the plane of the SET back- 
ground charge go and the amplitude Vin of incoming rf wave 
for g = 50, T = 0.01 eVCs, Ri:/Ro = 2000, and uj = wq. 
Vq = for panels (a) and (b), and Vo = O-Sc^/Ce for panels 
(c) and (d). 

tank circuit frequency = 1/ ^/LtCt, and the Q-iacXor 
Q = ^JLtJCt jR^- The operating point parameters are 
the efFcctive SET background charge go, dc bias voltage 
Vb, ampUtude Vi„ of the incident wave, and its frequency 
cj (in most cases we assume w = wo)- 

We use the SET parameters for normahzation, so that 
a natural unit for temperature is e^/Cs (where Cs = 
C\ +C2), the voltage unit is c/Ct., the RF-SET response 
dX/dqo (or dY/dq^) can be measured in units l/C^, and 
the unit for sensitivity dq^/s/KJ is e{RY,C-sY''^ (where 
i?s = i?i + i?2)- Notice that all considered magnitudes 
have a simple scaling with C^; however, there is no simple 
scaling with because of the dimensionless parameter 
Rh/Rq- 

Operation point optimization 

Figure 2(a) shows the RF-SET response dX/dqo on 
the plane of operating point parameters go and Vin for 
the case Q = 50, Rs/Ro = 2000 (i.e. i?s = lOOkfl), 
Vo = 0, Lu = ujo, and T = O.OleVCs (the SET is sym- 
metric, Ci = C2, Ri = R2)- Figure 2(c) is similar, except 



Vo = 0.5e^/Cs. We do not show dY/dqo because it is 
practically vanishing. Because of the same reason, the 
sensitivity 5qo/\/Af shown in Figs. 2(b) and 2(d) is cal- 
culated using only quadrature X. 

Both the response and sensitivity are obviously poor 
when the amplitude of oscillations at the SET is be- 
low the Coulomb blockade threshold Vt [this happens at 
Vin < Vt/2Q where V = (e/Cs)(l - 2go/e), and corre- 
sponds to the triangles in the lower left corners of Figs. 
2(a)-(d)], and they are also poor when Vin is much larger 
than this condition (notice that better response corre- 
sponds to larger dX/dqa, while better sensitivity corre- 
sponds to smaller Sqa). Even though the regions of rela- 
tively good response and sensitivity are similar, the max- 
imum response and optimum sensitivity are achieved at 
quite different points in the Vin — Qo plane. In particular, 
the amplitude Vin of the incoming wave is significantly 
larger for max;imum response, than for best sensitivity. 

In the present-day experiments, maximization of re- 
sponse [Figs. 2(a) and 2(c)] is still of the major impor- 
tance, because the noise from the next amplifying stage is 
still significant [if a large constant noise had been added 
to Sx in Eq. (42), then the sensitivity would be mainly 
determined by the denominator dX/dqo\. However, if 
the amplifier noise is small compared to the contribution 
from the SET shot noise, then the best operation point 
should optimize the shot-noise-limited sensitivity [Figs. 
2(b) and 2(d)]. In the case of comparable contributions 
from two noises, we have a trade-off between two regimes. 
In this paper we will concentrate on the analysis of the 
maximum response (MR) mode and optimized sensitivity 
(OS) mode, keeping in mind that experimentally optimal 
regime is somewhere in between, depending on the am- 
plifier noise. 

Figure 3 shows dependence of the response dX/dqo and 
sensitivity Sqo/^/Af on the dc bias voltage Vq in the 
MR and OS modes (optimizations are over Vin and go)- 
Several curves on each plot are for different Q-factors: 
Q = 10, 30, 50, 70, and 90, while other parameters are 
similar to the parameters of Fig. 2. One can see that the 
best response in both MR and OS regimes as well as the 
best sensitivity in OS mode are achieved at Vo = 0. This 
is because both positive and negative branches of the 
symmetric SET I — V curve (see Fig. 10 in Appendix) 
contribute equally at Vq = 0, and therefore the signal is 
maximal. 

In the MR mode at Vq = the optimum background 
charge [see Fig. 2(a)] is about go « 0.15 e (so the Coulomb 
blockade threshold Vj is about 0.7e/Cs) while the op- 
timum amplitude At of the SET bias VJ, oscillations is 
about l.le/Cs (these numbers have only weak depen- 
dence on Q, while the optimum Vin obviously depends on 
Q quite significantly). When Vq starts to increase, it be- 
comes advantageous to increase go (so Vt decreases) while 
Ah stays approximately constant, so that both positive 
and negative branches of the SET I—V still contribute to 
the response. However, since these branches cannot both 
contribute in the optimal way, the response decreases 
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FIG. 3. Dependence of (a) (b) the RF-SET response 
dX/dqo and (c)-(d) sensitivity Sqoj^J'Kf on the do bias volt- 
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qo are optimized for either maximum response [MR mode, 
panels (a) and (c)] or optimal sensitivity [OS mode, panels 
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Lj = ojo- Notice inoiiotoiiic worsening of the sensitivity with 
Q. The horizontal dotted lines in (b) show the low-Q low-T 
results^" corresponding to Eqs. (43) and (44). 



Dependence on Q -factor 



with Vq [see Fig. 3(a)]. For large enough Vq it becomes 
preferable to use only one (positive) branch [this regime 
corresponds to the lower maximum in Fig. 2(c), while the 
upper maximum corresponds to the two-branch regime] ; 
then the optimal wave amplitude drops to At ~ 0.2 e/Cs 
and the optimal qo corresponds to Vt slightly above Vq 
(by about 0.03 e/Cs). This change causes the kinks on 
the response curves in Fig. 3(a) and jumps down on the 
sensitivity curves in Fig. 3(c). Notice that in the MR 
mode the sensitivity at Vq — O.Se/Cs is better than at 
Vo=0. 

In the OS regime the optimal amplitude is significantly 
less than in the MR regime. It depends mainly on tem- 
perature, and for parameters of Fig. 3 the amplitude 
at the SET is typically between 0.08 and 0.1 in units of 
e/Cs- The "above blockade" voltage Vq + Ai, — Vt is few 
times smaller than Ai, and is comparable to the temper- 
ature. The best sensitivity is achieved at Vb = when 
both branches participate, while at large enough Vq when 
it becomes preferable to use only one branch, both sensi- 
tivity and response practically do not depend on Vq [see 
Figs. 3(b) and 3(d)]. 



The dependence on Q-factor is summarized in Fig. 4, 
which shows the response and sensitivity in the MR and 
OS modes for T = 0.01 eVCs, Rn/Ro = 2000, and 
ui = LOo- The results are presented for Vb = (thick 
lines), which provides the best MR and OS response and 
best OS sensitivity, and also for Vq = 0.5e/Cs (thin 
lines), which is a typical value for the case when only one 
branch of the SET I — V is involved. One can see that 
the response dX/dqo grows linearly with Q at small Q 
[see Eq. (11)] and reaches the maximum at Q around 50 
(this number is somewhat different for different regimes; 
for example it is almost 70 for the thin solid line), which 
is close to the crude theoretical estimate -\/RsA^ — 45 
for the impedance matching. However, unlike in the lin- 
ear model, this maximum does not correspond to the ex- 
act impedance matching. For example, the impedance 
matching (minimum of reflection) occurs at Q ~ 100 
for the upper curve in Fig. 4(a) and at Q ~ 80 for the 
curve second from the top, while for two lower curves 
(OS mode) it does not occur at all in a reasonable range 
of Q. 

In contrast to the response behavior, the RF-SET sen- 
sitivity monotonically worsens with increase of Q. Qual- 
itatively, this happens because the noise Sx has scal- 
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ing [see Eq. (35)], while the response dX/dqo has slower 
than Q dependence [see Eq. (11)]. This simple analysis 
predicts the parabolic dependence Sq^^ oc /{Rd/Ro), 
which crudely fits the curves in Fig. 4(b) using Rd/Ro 
about 1.5-3 times larger than R^/Rq, though the curves 
in Fig. 4(b) actually have slower than parabolic depen- 
dences at large Q. 

Comparing the MR and OS modes at Q = 30 (just 
some typical number) and Vb = 0, wc sec that the MR 
regime provides about 40% larger response, while the 
OS regime provides about twice better sensitivity. Even 
though these numbers depend significantly on the tem- 
perature and also depend on Q and Rs/Roy they show 
that the results for the MR and OS regimes are not too 
much different (not by an order of magnitude). 

Comparison of the cases Vq = and Vq = 0.5e/Cs 
shows that for the MR and OS response as well as for 
the OS sensitivity there is no much difference between 
these two cases, and the relative difference decreases with 
Q. In contrast, the difference between MR sensitivity at 
Vq = and at Vb = 0.5 c/Ce grows with Q and can 
become significant. 

The low-Q limit of the OS sensitivity for Vb = is well 
described by the formula^° (which also assumes T <C 
eV^E) 



0.020 



^go/^A/ ~ 2.65e(i?sCs)i/'(TCs/e')i/' (43) 

[we do not use a shorter formula 2.65 Cs{RsT)^^^ to 
emphasize natural normalizations], and similar limit for 
Vo = O.Se/Cs (one-branch case) is close to^° 

^go/^/A7 ~ 3.34 e{Rj:C^y^\TC^/e'')'/^ (44) 

[see Fig. 4(b)]. However, the theory of Ref.^" which as- 
sumes Q <C \/Rj:/Ro, is not able to describe significant 
change of the sensitivity with Q in Fig. 4(b). One can 
see that this dependence is even more significant in the 
MR mode. 



Dependence on temperature and SET resistance 

The numerical results on temperature dependence are 
shown in Fig. 5 for Ry:/Ro = 2000 and uj = ujq. Similar 
to Fig. 3 we show on four panels the response dX/dqo 
and corresponding sensitivity in the MR and OS modes 
(the optimal dc bias value Vb = is used). It is im- 
portant to notice that in the MR mode both response 
and sensitivity almost do not depend on temperature at 
T < 0.03 eVCs, and the RF-SET performance is still 
OK at temperatures ~ 0.1 e^/Cs (response and sensitiv- 
ity change less than twice compared to zero-temperat\irc 
case) . In the OS regime the response also has a very weak 
temperature dependence at T < 0.05e^/CE; however, the 
sensitivity is strongly temperature-dependent. The low- 
Q OS sensitivity can be accurately described by Eq. (43) 
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FIG. 5. Dependence of (a)-(b) RF-SET response and 
(c) (d) sensitivity on temperature T in the MR and OS 
regimes for several Q-factors. Vb = 0, R^/Ro = 2000, and 
u) = ujQ. Lowest dotted line in (d) is Eq. (43). 



up to temperatures ^ 0.05e^/Cs, [sec the lowest dot- 
ted line in Fig. 5(d)], and the curves for large Q-factors 
also follow the scaling Sqo/y/Af oc T^/^ at temperatures 
T < 0.05 e^/Cs. In the orthodox theory this T^/^ depen- 
dence is valid even at very small temperatures leading 
to infinitely good sensitivity {6qo/y/ Af 0); however, 
in reality the neglected contribution of cotunneling pro- 
cesses becomes significant in the OS regime at small T 
that changes the formalism.^^'^^ Comparing the results 
for MR and OS modes, we notice that while they are 
significantly different at low temperatures, the difference 
decreases with temperature, so that at T ^ O.le^/Cs 
the MR and OS results are already quite similar. 

Figure 6 shows the response and sensitivity depen- 
dence on the total SET junction resistance i?s for T = 
O.Ole^/Cs, Vb = 0, and uj = ojq (we show results only 
for R^/Rq > 500, because at i?E < 25kfi the theory is 
too inaccurate due to neglected cotunneling processes). 
In the case'^° Q <C ^/Rs/Ro the response dX/ dqo scales 
as R^^ [see Eq. (23)] and the sensitivity 6qo/y/KJ scales 

1/2 

as [see Eq. 43]. Correspondingly, the solid lines 

{Q — 10) in Figs. 6(a) and (b) are practically straight 
lines passing through the origin, and in Figs. 6(c) and 
(d) the solid lines are practically horizontal [the level de- 
termined by Eq. (43) is shown in Fig. 6(d) by the lowest 
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FIG. 6. (a)-(b) Inverse response and (c)-(d) sensitivity in 
the MR and OS regimes as functions of the SET resistance 
Rs for several Q-factors at T = 0.01, Vo = 0, and uj = ujo- 



dotted line, which fits well the solid line]. The depen- 
dence on i?E becomes nontrivial when Rs/Rq is compa- 
rable to Q^. It is interesting to notice that the inverse 

response remains a practically linear function of R-^ even 
for large Q, as seen in Figs. 6(a) and (b) (the slope of 
the lines decreases with Q while the offset increases). In 
particular, this means that in the orthodox model the 
RF-SET response continues to increase with the SET re- 
sistance decrease even when the matching condition is 
overshot. Decrease of Rs close or beyond the match- 
ing condition, worsens the sensitivity in comparison with 

scaling Sqo/\/Af oc R^^ [see Figs. 6(c) and (d)]. In the 
MR mode the sensitivity worsens with decrease of Ry: 
even in absolute units (e/VHz) when the matching con- 
dition is sufficiently overshot, while in the OS mode the 
sensitivity still improves with decrease of R^ in absolute 

1/2 

units, in spite of worsening compared to i?^ scaling. 



Effect of asymmetric biasing (Ci C2) 

Now let us discuss the effect of asymmetric effective 
capacitances of the SET junctions, Ci 7^ C2. Even when 
the physical junction capacitances arc equal, Cij = C2j, 
the effective capacitances Ci and C2 can be significantly 
different when the gate capacitance Cg is comparable to 
junction capacitances. In our model shown in Fig. 1 this 
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FIG. 7. (a)-(b) Dependence of the RF-SET response and 
(c)-(d) sensitivity on the dc bias voltage Vo for several 
values of the effective asymetry of the SET capacitances. 
T = O.OleVCB, Q = 50, J?s/i?o = 2000, and uj ^ loq- No- 
tice that the SET asymmetry ( "asymmetric biasing" ) slightly 
improves the RF-SET performance. 

happens due to asymmetry of signal source capacitances, 
Csi ¥= Cs2 [see Eqs. (1) and (2)]; in RF-SET experiments 
this effect is called asymmetric rf biasing of the SET.-'^-'^ 
For the conventional SET setup, the biasing asymmetry 
is not important at all, because of the formal SET equiv- 
alence under transformation^^ 

Ci^Ci + AC, C2^C2- AC, go ^ go - VbAC (45) 

for arbitrary AC, that means that the asymmetric bi- 
asing (different Ci and C2) can be simply corrected by 
the background charge shift. However, for the RF-SET 
there is no simple equivalence because the bias voltage Vb 
changes in time. (Effective capacitance asymmetry can 
still be controlled by addition of the extra rf signal to 
the SET gate.) Therefore it is interesting to find if the 
asymmetric biasing (Ci ^ C2) is better or worse than 
the symmetric case (Ci = C2). 

Figure 7 shows the Vq dependence of the RF-SET re- 
sponse and sensitivity in the MR and OS modes for sev- 
eral ratios C1/C2 (with fixed total capacitance Ci + C2). 
One can see that in the asymmetric cases the best MR 
response and OS sensitivity are achieved at nonzero Vq 
and are typically better than the corresponding values 
for the symmetric case. Therefore, the asymmetric bias- 
ing is actually preferable (that eases the concern about 
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this issue discussed in Ref.-'^'^); however, the advantage 
is rather minor. In particular, the OS sensitivity in the 
asymmetric case is still limited by Eq. (44) ; therefore the 
possible improvement is less than the sensitivity decline 
due to high Q-factor [see Fig. 4(b)], and even less because 
of the difference between Eqs. (43) and (44). 



Carrier frequency detuning from the resonance 

So far we have considered the resonant case lo = ujq. 
[Actually, the exact resonance is at the frequency lo — 
a;o(l — l/2Q|^y) which is very close to wq-] In this case 
the quadrature component X is much larger than Y (at 
Q ^ 1) and therefore the RF-SET response and sensi- 
tivity in respect to monitoring X quadrature practically 
coincide with that for monitoring the reflected wave am- 
plitude A (we denote with A the amplitude Ai of the 
first harmonic) or monitoring the optimized phase-shifted 
combination X* = Xcos(fi + Ysuv(f. The detuning of 
(jj from the resonant frequency ljq leads to a significant 
magnitude of Y quadrature, and so to different results 
for different ways of reflected wave monitoring. 

For a small frequency detuning Aw = w — ujq the 
simple linear analysis using Eq. (11) predicts the non- 
monotonic shape of the X-rcsponse frequency depen- 
dence with zeros at w/wq = 1 ± 1/2(5l and with the 
full width at half height (FWHH) equal to O.SSwq/Ql; 
while for the amplitude monitoring it predicts FWHH of 
VSwo/Ql. [Eq. (8) which assumes both Rd and Vin being 
frequency- independent, gives the response FWHH equal 
to (v^— 2)^/^a;o/Qz, and uq/Ql ior X and A monitoring, 
respectively]. In reality the situation is more complex 
because even neglecting the nonlinear contributions, the 
effective SET resistance Rd (and therefore Ql) depends 
on Vb{t) amplitude, which depends itself on w, Rd, and 
Vin- Moreover, the optimized point in the Vin — Qo plane 
is also frequency-dependent. 

Figure 8 shows the numerical results for the frequency 
dependence (in the vicinity of the resonance) of the op- 
timized response and sensitivity. Each panel shows four 
curves, which correspond to monitoring the components 
X (thick solid line), Y (dashed), optimized phase-shifted 
combination X* (thin solid), and amplitude A (dotted 
line). At each frequency we perform optimizations of 
the response (MR mode) and sensitivity (OS mode) over 
Vin — Qo plane for all four monitored magnitudes, so the 
optimizations are different for different curves. The Q- 
factor is equal to 50, which is close to the estimate of the 
impedance-matching value (see also Fig. 4); however, in- 
stead of the naively expected value Ql ~ (5/2 = 25 for 
the loaded Q-factor, it is = 40 in the MR mode since 
QsET — 199 (these values are calculated a.t uj = u>o using 
effective SET resistance). The shape of the X-response 
dependence in panel (a) is close to the prediction of Eq. 
(II) using Ql = 40, though the minima do not reach 
zero and the curve beyond the minima is shifted up. The 



w 

1/5 
PL, 



0.025 - 


MR mode 


(a): 


0.020 - 


Vo=0 




0.015 - 


r \ I' 

/ */ Ir 


\ ■ 

\\ 


0.010- 


/ |i '1 


\ \ 


0.005 - 


\ / 1 ' 


/ ^ 
1 \ ~ 


0.000 







OS mode (b) : 

\dYldq,\C^ - 

dX Idq^ ; 

\dA/dqo\C^ - 



0.95 



1.00 

Co/Wf, 



1.05 



O 0-8 

O 

0.2 







(c) 










MR mode " 



0.95 



1.00 

m/con 




1.05 



FIG. 8. Frequency dependence of the (a)-(b) RF-SET re- 
sponse and (c)-(d) sensitivity in the MR and OS modes. 
Each panel shows four curves corresponding to monitoring 
of X quadrature (thick solid line), Y quadrature (dashed), 
their optimized phase-shifted combination X* (thin solid), 
and the first harmonic amplitude A (dotted). Q = 50, 
T = O.OleVCs, Rs/Ro = 2000, Vb = 0. 



MR y-response is practically zero at uj = ujo (since Y 
vanishes); however, it becomes comparable to the reso- 
nant X-response at some frequency detuning; the overall 
shape is close to the prediction of Eq. (11), but the max- 
ima are about 15% higher. The MR A-response curve 
is about 50% wider than the prediction of the linear 
theory (using Ql = 40) and is very close to the curve 
for monitoring optimal phase-shifted quadrature X* . As 
expected, the A- and X*-responses are better than X- 
and y-responses at finite detuning (for both MR and OS 
modes). 

The sensitivity (both OS and MR) for X* or A mon- 
itoring at finite detuning is also better then for X or 
Y monitoring. An interesting observation is that while 
the RF-SET response decreases with detuning, the sen- 
sitivity slightly improves with detuning for X* and A 
monitoring. This effect is similar to the sensitivity im- 
provement with the decrease of Q-factor (see Fig. 4). 
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wave is in resonance with the tank circuit. T = 0.01 e^/Cs, 
Rj:/Ro = 2000; Vo = O-Sc/Ce for uj = u>o/2 and Vo = for 
ui = oJo/3. 



Monitoring of resonant overtone 

In experiments the incoming rf wave is usually tuned 
close to the resonance with the tank circuit; in this case 
the contribution of overtones in the reflected rf wave is 
small in spite of significant SET I — V nonlinearity (the 
SET nonlinearity has been recently used^^ for rf mixing) . 
However, if the nth overtone is in resonance, uj w ujo/n, 
then the reflected wave may have a significant contribu- 
tion from this overtone, and the RF-SET operation can 
be based on monitoring this overtone. ''''' The use of differ- 
ent frequencies for the incident and reflected waves may 
be advantageous for some applications. Also, it may be 
useful to have the absence of the monitored reflected wave 
when the SET is off (no current), while in the conven- 
tional regime this case corresponds to the largest reflected 
power. One more possible advantage is somewhat easier 
control of the amplitude of the SET bias voltage oscilla- 
tions, since now it is more directly related to the incident 
amplitude Vin, while in the usual regime the relation de- 
pends on Ql [see Eq. (27)] which varies with operating 
point. (The disadvantage is that the incident amplitude 
Vin should be much larger than in usual regime that may 
lead to heating problems.) 

Figure 9 shows the RF-SET response and sensitivity 
for (jj = uJo/2 and w = wq/S, in respect to monitoring I2 
and I3, correspondingly (the X-quadratures are small). 



We use Vo = in the case iv = wo/3 and Vq = 0.5e/Cs in 
the case u = luo/2 [for Vo = there is no second overtone 
because of the I — V curve symmetry even at nonzero go 
- see Fig. 10(a) and Eq. (17)]. As seen from Fig. 9, the 
MR responses and OS sensitivities in the two cases (sec- 
ond and third resonant overtone) are quite similar, that 
is related to a strong nonlinearity of the SET I—V curve. 
Moreover, comparing Figs. 4 and 9 (the parameters are 
the same) we sec that the RF-SET performance in the 
regime of a resonant overtone is comparable to the per- 
formance in the conventional regime ui = uiq (the MR 
response and OS sensitivity are worse by only about 1.5 
times). Combined with the advantages discussed above, 
this makes the regime of resonant overtone potentially 
useful in experiments. 



IV. CONCLUSION 

In this paper we have developed the formalism for the 
calculation of the response and shot-noise-limited sensi- 
tivity of the normal-metal RF-SET, extending the theory 
of Ref to the case of arbitrary large Q-factor of the tank 
circuit and arbitrary frequency of the incident rf wave. 
The formalism has been used to analyze numerically the 
dependence of the RF-SET response and sensitivity on 
the operation parameters. 

We have mainly studied two operation modes, opti- 
mized over the rf wave amplitude and the SET back- 
ground charge, corresponding to either maximum re- 
sponse (MR mode) or optimum sensitivity (OS mode). 
The rf amplitude for the optimum sensitivity is typi- 
cally significantly smaller than for maximum response. 
The MR mode is the best experimental regime when the 
preamplifier noise is more significant than the effect of 
the SET noise, while the OS mode is the best when the 
preamplifier noise is negligible. Analyzing the perfor- 
mance dependence on the SET dc bias voltage, we have 
found that the best response and sensitivity are achieved 
at zero dc bias, though finite biasing docs not change the 
RF-SET performance much as long as it is within the 
Coulomb blockade range. 

We have found that the dependence of the RF-SET re- 
sponse on the unloaded Q-factor of the tank circuit has 
a maximum at Q comparable to the simple impedance 
matching estimate y/Rs/Ro. In contrast, the RF-SET 
sensitivity monotonically worsens with increase of Q; the 
dependence is approximately parabolic and the sensitiv- 
ity change can reach few times (compared to the low-Q 
case) if Q is chosen too large. This means that to im- 
prove the sensitivity in an experiment, it is better to 
"undershoot" the Q-factor compared to the impedance 
matching case, than to "overshoot" it. 

Studying the temperature dependence we have found 
that the RF-SET response saturates approximately at 
temperatures T < 0.03e^/Cs, that translates into 200 
mK for a typical capacitance value Cs = 300 aF. The op- 
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timized sensitivity continues to improve as T^/^ at lower 
temperatures until it reaches the quantum limitation due 
to cotunneling^^'^^ (not studied here). The "orthodox" 
RF-SET performance improves with decrease of the SET 
resistance; however, for high Q-factors the dependence is 
significantly slower than the low-Q scaling R^^ for the 

response and R^"^ for the sensitivity. 

We have analyzed the effect of the asymmetric rf bias- 
ing of the SET which leads to unequal effective capaci- 
tances of the SET and found that such asymmetry does 
not worsen the RF-SET performance (even slightly im- 
proves it). This answers the concern about asymmetric 
rf biasing raised in Ref.^^. 

We have analyzed the effect of the carrier frequency de- 
tuning from the resonance and found that the decrease 
of the RF-SET response with detuning can crudely be 
described by the simple linear theory (though difference 
in linewidth can reach 50%); however, the estimate of the 
loaded Q-factor determining the linewidth is not simple 
since it significantly depends of the operation point. Even 
with the frequency detuning, the RF-SET performance 
for monitoring the reflected wave amplitude (by rectifica- 
tion) is found to be very similar to the case of optimal ho- 
modyne detection (one-channel mixing with the optimal 
phase). It is important that the mixing phase optimizing 
the sensitivity can be different from the phase optimizing 
the response. Unexpectedly, in contrast to the response 
decrease with detuning, the optimized sensitivity slightly 
improves with detuning. 

We have also analyzed the operation regime for which 
an overtone of the incident rf wave is in resonance with 
the tank circuit, and found that the RF-SET perfor- 
mance in this regime is comparable to the performance 
in the conventional regime. Taking into account an ad- 
vantage of the frequency separation between the incident 
and monitored waves, this operation mode may be ex- 
perimentally useful. The theoretical proposal of the res- 
onant overtone regime"^^ has been recently realized exper- 
imentally in the group of Keith Schwab. Experimental 
RF-SET sensitivities using the second and third resonant 
overtones have been found practically coinciding with the 
sensitivity in the conventional regime. 

Now let us compare our theoretical results for the RF- 
SET sensitivity with the experimental value of 9 ^e/VUz 
from Ref.^^ for the normal- metal case (the sensitivity 
of the superconducting RF-SET was significantly bet- 
ter: 3.2 /ie/\/Hz). Using the experimental parameters 
Ce = 267 aF, i?E = 43 kil, assuming the temperature 
T = 70 mK, and using Fig. 6, we get the MR sensi- 
tivity of 3.2/ie/v^ for Q = 30 and 5.8^e/\/m for 
Q = 50, while the OS sensitivities are 1.3/ie/\/Hz for 
Q = 30 and 1.9 ^e/VHz for Q = 50 [notice that all these 
numbers are significantly higher than the low-Q OS es- 
timate 0.9/ie/\/Hz using Eq. (43)]. Unfortunately, it is 
not possible to extract the unloaded Q-factor value from 
Ref.^^; however, it seems to be between 30 and 50. It 
is also not known if the RF-SET operation point was 



closer to the MR or OS mode (we guess the response was 
more important). Nevertheless, we see that the differ- 
ence between the theory and experiment is few times. We 
guess that the difference is mainly due to the preampli- 
fier noise. Some contribution may also come from higher 
effective tcmperatiiro than we assumed and from a non- 
optimal operation point. Another contribution to the 
difference may come from the neglected here effect of co- 
tunneling, which limits the sensitivity. However, a rough 
estimate^^ of this limit 5qo/^Af ~ y/hC^ gives the value 
of 1 /Lte/\/Hz, so it is unlikely to be the major reason for 
the difference. Wc hope that the further experimental 
progress will bring the RF-SET sensitivity really close to 
the theoretical limit. 

There are still many theoretical questions about the 
RF-SET performance, not answered in this paper. For 
example, as seen from the above estimate, the account of 
cotunneling and study of the quantum operation of the 
RF-SET is starting to be important for experiments. The 
development of the theory for superconducting RF-SET 
is even more important since in the majority of experi- 
ments with RF-SETs the superconducting state is used. 
There is still no rigorous theory of the frequency depen- 
dence of the RF-SET sensitivity. It is also important to 
consider the backaction from the RF-SET and analyze if 
it can in principle be used as a quantum detector with 
high quantum efficiency (ideality). These problems are 
going to be the topics of further studies. 
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APPENDIX: CURRENT AND NOISE 
CALCULATIONS FOR THE SET 

The SET current / and its low frequency shot noise Si 
have been calculated using the orthodox theory of single- 
electron tunnehng,^ assuming that a stationary state of 
the SET is achieved at any moment of time (i.e. //e <C 
ui). We assume normal metal (non-superconducting) case 
and calculate the tunneling rates (™) of electron tun- 
neling to {+) or from (— ) the island through the first or 
second junction as^ 



(m) 



e^Rj[l - exp(-VF/(m)/r)] 



e 



,.jVbC^ 1 go 



(Al) 
(A2) 



where m is the number of extra electrons on the SET, 
j = 1,2 denotes the junction, and T is the temperature. 
The current [Fig. 10(a)] is calculated as 
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I = eJ^[T+{m)-T^{m)]ast{m), (A3) 

m 



where the stationary probabihty distribution of the 
charge states <Tst {m) satisfies equations 

ast (m) ^ r+ (m) = a.t (m + 1) ^ T" (m + 1) (A4) 
J j 

and J2^<ystim) = 1. 

The low-frcqucncy spectral density Sj of the SET cur- 
rent [Fig. 10(b)] can be calculated as^^ 

Si = -Ae^ \^t{m') - T^{m') - 7/e](f-i)„,„ 

X [r+(TO - l)ast(.m -l)-Tr^{m+ l)G,i{m + 1) 
-(//e)a,*(m)] +2e2^a,t(m)[r+(m) +rr(m)], (A5) 

m 

where T is the three-diagonal matrix of the SET charge 

evolution, Tkm = I]j[^m,fc-ir+(m) + 8m,k+i^~ [m) - 
<5rri,fe[r^(™) +rj(m)]]. Notice that Eq. (A5) can be eas- 
ily used numerically even though the matrix T is singular 
and therefore does not have a unique inverse; Eq. (A5) is 
constructed in a way that the non-uniqueness is not im- 
portant, and therefore the standard algorithm for solving 
a linear system of equations with three-diagonal matrix 
can be readily used. (There is no problem of the matrix 
singularity at a finite frequency Wg, when the inverse of 
the matrix F — iwsl should be calculated.) Actually, in- 
stead of using Eq. (A5), we have used a somewhat faster 
algorithm for calculation of Si, described in Section VII 
of Ref.23. 
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